General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 





DEPARTMENT OF ELECTRICAL ENGINEERING 
SCHOOL OF ENGINEERING 
OLD DOMINION UNIVERSITY 
NORFOLK, VIRGIN U 



AN INVESTIGATION TO MODEL RADAR BACKSCATTER 
FROM RIDGED FIRST YEAR SEA ICE 



UJ 

CO 

UJ 


CV 

> 



By 

Betty Lou W. Jackson 
and 

John W. Stoughton, Principal Investigator 


Final Report 

For the period ending December 31, 1981 





Prepared for the 

National Aeronautics and Space Administration 
Langley '’^search Center 
Hampton, irginia 


I 


Under 

Research Grant NASl-15648 
Task Authorization No. 51 
William L. Grantham, Technical Monitor 


(NA3 A- 

Ck- I 730 5 0) AN INVE 

f ACAu 

E A CK o V, / A T E F lA V i ■ j- 

ICE FO 

I- Hit P FH IvO EN L'l N J 

Fiiid 1 

e port (0 i(! P' i a 

Vd. ) 

9 3 p iir 05/Mr AC' 1 


lI.iAlICN :c 
( I; L H r n" Y lit I ^ r A 
L iH; K ? r L 1' 3 1, 13 3 1 

U:.i v<, , ^ oik , 

i;3CL 20 N G3/3 2 




oci 


\ f'' 




access 


0€Pt. 






^ .3 3 - 3 0 1 9 7 
U I. c i a t 


August 1983 



DEPARTMENT OF ELECTRICAL ENGINEERING 
SCHOOL OF ENGINEERING 
OLD DOMINION UNIVERSITY 
NORFOLK, VIRGINIA 


AN INVESTIGATION TO MODEL RAD^ BACkSCATTER 
FROM RIDGED FIRST YEAR SEA ICE 


By 

Betty Lou W. Jackson 
and 

John W. Stoughton, Principal Investigator 


Final Report 

For the period ending December 31, 1981 


Prepared for the 

National Aeronautics and Space Administration 
Langley Research Center 
Hampton, Virginia 


Under 

Research Grant NASl-15648 
Task Ai’.chorization No. 51 
William L. Grantham, Technical Monitor 




Submitted by the 

Old Dominion University Research Foundation 

P.O. Box 6369 

Norfolk, Virginia 23508 


August 1983 


TABLE OF CONTENTS 


PAGE 

LIST OF TABLES iv 

LIST OF FIGURES v . 

SYMBOL DEFINITIONS vli 

SUMMARY viii 

Chapter 

1» Introduction 1 

1.1 Forward 1 

1.2 Need for Ice Research 1 

1.3 Need for Data About Ice Ridges 3 

1.4 Measureaent Technlqiies 4 

1.5 Research Objective 7 

2. Instrumentation and Data Processing 9 

2.1 Introduction 9 

2.2 Mapping Camera li 

2.3 Laser Profiler 16 

2.4 Microwave Scatterometer 18 

2.5 Data Processing 20 

2.5.1 Instrument Perspective Differences .... 20 

2.5.2 Data Selection 21 

2.5.3 Surface Slope Processing 22 

3. Backseat ter Model 31 

3.1 Introduction 31 

3.2 Backscatter Signatures 31 

3.3 Backscatter Model 38 

3.3.1 Calculation of Incidence Angles 41 

3.3.2 Instantaneous Reflected Power Model ... 43 

3.3.3 Average Reflected Power Model 44 

3. 3. 3*1 Equivalent Resolution Cell ... 45 

3. 3. 3. 2 Antenna Weighting 48 

3.3.3. 3 Calculation of Area 51 

3.3.4 Addition of Radar Measurement Noise ... 51 

4. Model Simulation 33 

4.1 Introduction 53 

4.2 Data Selection 54 

4.3 Adjustment to • 54 


iii 


PAGE 


4.4 Time Adjustment by Features • 56 

4.5 Adjustment to Signature Slopes 57 

4.6 Stereo Analysis 61 

5. Conclusions 74 

5.1 Concluding Remarks 74 

5.2 Future Research 75 

References 78 


^pendlx A 


81 



Iv 


LIST OF TABLES 

W Stereo Profile Statistics 


PAGE 

70 


V 


LIST OF FIGURES 


FIGURE page 

2.1 Flights Lines for SIRE 79 .....«••• • 10 

2.2 Relative Location of Instruments' Illumination Areas . 12 

2.3 Facsimile of Aerial Photogram with Young Ice ... 13 

2.4 Facsimile of Aerial Photogram of Multi-Year Ice ... 14 

2.5 Laser Profiler 17 

2.6 Scatterometer Antenna 19 

2.7 Physical Representation of Phase Shift 24 

2.8 Removal of Phase Shift from Data • 25 

2.9 Aircraft Motion In Data ........... 27 

2.10 High-Pass Filter Attempt to Remove Aircraft Motion . . 28 

2.11 Aircraft Motion Removal 79 

3.1 Onstott's Signature Curves 32 

3.2 Illustration of the Relationship Between 

0 and 6jj for Flat Ice 34 

3.3 Illustration of Equal Return for ©n ^ ..... 35 

3.4 Facsimiles of Surface Photographs of Ridges .... 36 

3.5 Facsimile of Aerial Photogram of Ridged Ice .... 37 

3.6 Illustration of Facets in Equivalent Resolution Cell . 39 

3.7 Flowchart of Modeling Processes ........ 40 

3.8 Illustration of(^“6u+a 42 

3.9 Flowchart of averaging Processes 46 


Vi 


3.10 Illustration of Overlapping Resolution Cells .... 47 


3.11 RADSCAT Antenna Pattern 49 

3.12 Illustration of Weighting Sequence 50 

4.1 Selection of Data from Photograms 55 


4.2 Illustration of Correction for Feature Alignment ... 58 

4.3 Plot of Initial Predicted Return and RADScAT Series . . 59 


4.4 Plot of Final Predicted Return and RADSCAT Series . . 62 

4.5 Illustration of Stereo Profile Location 

Within Resolution Cell 64 

4.6 Facsimile of Aerial Photogram for Region A 65 

4.7 Facsimile of Aerial Photogram for Region B 66 

4.8 Plots of Stereo Profiles for Region A 67 

4.9 Plots of Stereo Profiles for Region B 68 


4.10 Plot of Decorrelation Time Between Profiles .... 69 


4.11 Plot of Predicted Return Based on 

Stereo Profiles for Region A 72 

4.12 Plot of Predicted Return Based on 

Stereo Profiles for Region B 73 

A.l Geometry of Scatterometer Footprint 82 


I 

I 



A. 2 Flowchart for RADSCAT 


83 


vil 


B 

br, bf 
d 

h(l) 

i 

j 

k 


Dr I ®f 


P(k) 

PR/S 

Pslm 

W(0) 

w(l) 

a(l) 

<)>k 

a** 

a[4»] 

OdB[4>I 

cr(()>it) 

6 

6a 


Symbol Definitions 

the 3dB beamvldth of the RADSCAT antenna 

Intercepts of the linear approximations of the signature 
curves for rough and flat Ice respectively 

distance between laser samples 

height as measured by the laser profiler 

Index of facets In the equivalent resolution cell* 100 Hz 

t 

Index of RADSCAT Integrated resolution cells, 10 Hz 
Index of angle bins 

slopes of the linear approximations of the signature 
curves for rough and flat Ice respectively 

percentage of the facets of the equivalent resolution 
cell In the k^h bln 

measured total return power In milliwatts, from the 
RADSCAT 

simulated total return power In milliwatts, the output 
of the model 

reflected power from surface due to antenna pattern and 
Integration time 

dlcrete weighting sequence used to approximate W(6) 
surface slope of the 1^^ facet 

true Incidence angle of the 1^^ facet, measured from 
the Incident wave to the surface normal 

the average angle of the k^^ angle bln 

normalized radar cross section, from Onstott's curves 

strength of return, proportional to rf , used to 
calculate pgim. In milliwatts 

10 logio C7[4>1 

average strength of return over the k^^ range of angles 

angle measured from boreslght of antenna 

the elevation angle of the antenna measured from nadir 


vlll 


AN INVESTIGATION TO MODEL RADAR BACwCATTER 
FROM RIDGED FIRST YEAR SEA ICE 

By 

Betty Lou W. Jackson^ and John W. Stoughton^ 

SUMMARY 

This paper investigates the response of a radar type sensor, the micro- 
wave scatterometer, to meter scale roughness of first-year sea ice. The 
scatterometer measures the absolute backscattered potfi'.r from target, and use 
of the radar equation permits calculation of the surface reflectivity 
expressed as the radar scattering coefficient, o*, also called the normal- 
ized radar cross section. 'Diis project develops an empirical, one-dimen- 
sional model to predict the back-scattered power from ridged first-year sea 
ice of known surface topography, by interpreting the surface as a series of 
angular facets and summing the proportional power from each facet. The 
model is evaluated by comparing the predicted values to scatterometer meas- 
surements obtained during the National Aeronautics and Space Administra- 
tion's (NASA's) Sea Ice Radar Experiment (SIRE) in 1978 and 1979. 

The input to the model is surface slope, which for SIRE was calculated 
from the ice topography, as measured by an airborne laser profiler. By 
limiting our consideration to first-year ice, the modeled radar backscatter 
is a function of slope only. The correlation between the predicted and the 
measured power is .62-. 68, and it is concluded that the model adequately 
describes the backscatter from regions of first year sea ice whose dimen- 
sions are of a few tens of meters. 


^Graduate Research Assistant, Department of Electrical Engineering, School 
of Engineering, Old Dominion University, Norfolk, Virginia 23508. 

^Associate Professor, Department of Electrical Engineering, School of 
Engineering, Old Dominion University, Norfolk, Virginia 23508. 


CHAPTER ONE 


INTRODUCTION 


1*1 Background 

Sea lce» the frozen surface of sea water, has a significant 

effect on our lives* At its maximua extent, this frozen layer covers 
+6 2 

40.6 X 10 km , nearly thirteen percent of the world's oceans. 

Sea ice plays a critical role in maintaining the earth's heat balance 
and directly Influences the climate* It is an obstacle to shipping 
and to the recovery of an estimated 200 billion barrels of oil and 300 
trillion cubic feet of natural gas [1]. Yet, despite its importance, 
less is known about sea ice than about any other part of the earth's 
surface. Only in this century have scientific ice expeditions reached 
the Arctic, and most of our present knowledge of sea ice has been 
acquired within the last two decades. To date, the properties of sea 
ice are not clearly described, much less understood. 

1*2 Need for Ice Research 


The study of sea ice has become interdisciplinary, fcom fore 
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casting vesthsr to mapping ocean currents and Ice drift* Weeks [2] 
points out that meteorologists are Interested In the patterns and 
locations of Ice ridges and leads (open cracks In the Ice sheet) 
because ridging and Ice break-up become most severe In the period 
Immediately preceding storms, and because ridges have a direct effect 
on the wind . Meteorologists suspect also that the location of the Ice 
edge Influences the frequency and paths of cyclones* Climatologists 
need data on surface temperatures, heat transfer and atmospheric 

pressure at high latitudes [3]* Oceanographers can study surface 
water currents by trading ice floes* As reported by Weeks [4], they 
need Information about, the dispersion of dense melt water and about 
the effects of heat and salinity on ocean currents* Biologists are 
studying sea Ice In an attempt to understand how the seasonal changes 
In Ice temperature and salinity affect the life forms of the arctic 
area, and whether the Ice biota affect the formation of the Ice [4]* 

In addition to the scientific uses, there are commercial and 
industrial applications for Ice Information as well* Shipping could 
be speeded end Its cost reduced If Ice modelers could predict regions 
where leads will occur which would facilitate Improved navigation 
within the Ice zone* The design of Ice breakers would be Improved If 
the shipbuilding Industry had realistic statistics for the average 
thickness and age of ice, and the number of ridges per nautical mile 
[5]* The petroleum Industry requires accurate tactical Ice fore- 
casting since an unexpected Ice movement or breakmp could demolish 
offshore oil rigs [6], and grounded keels (underwater protrusions of 
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ridg«$) of laod-fAtc rld^oc could fcour the shore end seebotcoB 

ceusiog extensive desMge to cebles end pipelines [7J. For defense 
purposes, the Nevy hes been studying the hydroacoustics of ice-covered 
water and the effects cf ice keeis on sooab [4J. Of late, there has 
been interest in icebergs as a source for fresh water for arid 
countries. It should be apparent that the need for ice research is 
real, and that it goes beyond oiaple scientific curiosity. 

1 . 3 Ne ed for Data About Ice Ridge.i 

The characteristic of sea ice which is of concern in this 
thesis is ice ridging, which is loosely defined as a linear distri- 
bution (Bound) of ice broken and buckled by coapreasive forces, and 
having heights of approxiaately 1/2 to several aeters. As aentloned 
previously, the shipping and shipbuilding Industries arc significantly 
affected by ridges: their height, width, keel depth, location, orien- 
tation and distributions being the aajor ridge paraaeters of 
laportance. The need exists for Inforaation on these paraaeters in 
real tlae and even future tlae (forecasting) for navigation, and over 
five to ten year tiaespans for ice breaker design. 

Another application for ridge data is ice dynaalcs aodellng. 
Ridges fora where leads have previously occurred; as a result, by 
utilizing past years' data, aodclers can test their predictions of 
where leads are likely to occur. Because the water in open leads 
peraits greater heat transfer than the surrounding ice, the location 
of leads is an iaportant input to cliaatology as well as to weather 
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forecasting. The wind is s factor in ti:* direction and velocity of 
ice Boveaent and is therefore a factor in ridge fonsation. The 
ridges, in turn, affect the vind drag of the ic«, changing the force 
which produced thea [2]. A better understanding of * he effect ridges 
and leads have on the wind would enable aeteorologists to make sore 
accurate weather forecasts, and the input of ridge location to ice 
forecasters' models would allow better prediction of ice novenent [8]. 

1.4 Measureaent Techniques 

Due to the inhospitable climate and the difficult access to the 
ice fields, the Inforaatlon available to researchers has been 
limited. Prior to the 1960's the ice data set consisted only of 
observations made from ships, sparse measurements made from surface 
based camps, and interpretations of aerial photograas (scaled 
photography). After considering some of the Ice Information required, 
such as surface temperatures, ice salinity and density. Ice age and 
thickness, lead patterns, and ridge distributions. It becomes apparent 
that these three aetlnods of collecting Ice data are not capable of 
producing a synoptic view of the Arctic. Ships are restricted to the 
open water; otherwise they nm the risk of being frozen Into the Ice 
mass. Reinforced vessels, which would be able to enter the i.ce field, 
^.re prohibitively expensive. Scientists working from surface based 
camps located farther into the Ice pack can make er.tenslve 
measurements from an area up to 10 In square. The surface data and 
lea cores they collect can provldn valuable, but very localized. 


information on ice structure. Through aerial reconnaissarce vlder 
areas can be covered, but since the Arctic is in constant thermal 
flux, the area is frequently cloud-covered and turbulent, which limits 
the opportunities for flight. The long periods of arctic daricness 
allow little aerial photogramme try even when the weather permits. An 
ice monitoring system needs to be developed which can be used over the 
expanses of the Arctic which is functional during poor weather and 
darkness, and which is preferably not surface based. Studies are 
being conducted on an alternate technique which can meet all of these 
requirements; that technique is remote sensing by satellites. 

Many remote sensors are under continued investigation to 

« 

determine their feasibility for ice monitoring and are being developed 
for satellite application based largely on encouraging results from 
aircraft experiments. Weeks summarized the current state of remote 
sensing [9]. Imaging, visible sensors, such as the Very High 
Resolution Radiometer (VHRR) on board the National Oceanic and 
Atmospheric Administration's (NOAA's) satellites, and the 
multlspectral scanning system flown on the Landsat satellite, have 
shown the capability to resolve large scale features (80 - 800 m) and 
thereby permit the tracking of ice floes and the Icepack. Lower 
resolution thermal sensors, also on board the NOAA satellite, are 
being used to define large leads and ice extent by differentiating 
between the temperature of ice and the surrounding water. Whereas the 
optical and infrared remote sensors require clear skies (and suitable 
lighting for the visible), the use of passive microwave sensors allows 


penecration of moderate cloud-cover as well as day and night operation 
to distinguish between ice and water as well as between ice types 
(ie., multi-year, first-year and thin ice). This technique uses 
sensitive microwave receivers (radiometers) to measure a gray-body 
electromagnetic radiation from ice and water which differ because of 
their emissivities and temperatures. Microwave radiometers have been 

deployed on Nimbus-5, 7 and 8 s&tellltes. Among active microwave 

• 

systems, there are studies on radar altimeters, synthetic aperture 
radar (SAR) and scatterometers. Radar altimeters, flown on Geos-3 and 
Seasat satellites, are nadir-looking radars which measure the height 
above the surface with a resolution of a few centimeters. For ice 
remote sensing, the altimeters use the difference in backseat tering 
amplitude to distinguish the ice/water boundaries. SAR and 
scatterometry both measure the strength of the radar return to provide 
roughness information, whether the surface is ice or water. SAR 
employs digital co. >lat <on processing to synthesize the large antenna 
necessary to provide a good resolution image (approximately 25 km) 
from space. Ice parameter extraction algoritnms for SAR are based on 
photo-interpretation techniques. The scatterometer is a quantitative 
active microwave sensor and has the ability to discriminate between 
water and ice and between multi-year and first-year ice over coarser 
resolution of tens of kilometers from satellite altitudes. Both SAR 
and scatterometer were on-board Seasat. 

Of these remote sensors, it is the data from the microwave 
scatterometer which will provide the comparison data set for the model 
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developed in this thesis* Further, the laser profiler and the mapping 
camera will be used as auxiliary sources of data to calculate the 
scatterometer response to ice roughness* 

1*5 Research Objectives 

Specifically, in this thesis, an empirical model is developed 
to predict the backscatter from ridged first-year sea ice, given the 
one-dimensional surface slope as the only input parameter* Experi- 
mental data are used where the slope is determined from the ice 
surface topography provided by the laser profiler, and comparisons are 
made between the calculated backscattered power and the value measured 
by the scatterometer to establish the validity of the model. The 
premise is, that if the above comparisons are highly correlated, then 
the model can be used to simulate the response of a satellite radar 
sensor, and thereby evaluate if the scatterometer can be used over 
first year sea ice to measure the degree of ridging* It is not within 
the scope of this thesis to reverse the modeling process and develop 
an algorithm to quantify the roughness based on scatterometer data. 

To this author's knowledge, investigations such as this have 
not been made* Models have been developed using existing data to 
describe and predict the physical phenomena of ice, such as ridge 
height or drift direction [10,11]* Further, theorists have been 
attempting to develop electromagnetic models to predict the volume and 
surface scattering of sea ice [12,13]. However, we know of no work 
where the response of a radar sensor to ice ridges has been modeled* 
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The Instruments used In this thesis are discussed In Chapter 
Two, along with the significant details of the Sea Ice Radar 
Experiment aircraft missions of 1978 and 1979 (SIRE)* Data processing 
for the laser is also discussed in that chapter* 

The modeling processes will be described in Chapter Three, and 
will include the development of the model* Results of the model 
simulation will be presented in Chapter Four* The predicted versus 
observed backscatter correlation will be determined, as well as the 
results of tests for statistical significance* Conclusions and future 
research implications will be given in Chapter Five* 
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CHAPTER TWO 


INSTRUMENTATION AND DATA PROCESSING 


2.1 Introduction 

In response to the need for remote sensing of Ice, as discussed 
In the Introduction, the National Aeronautics and Space Administration 
(NASA) has conducted a number of microwave remote sensing 
experiments. This thesis uses data acquired on two aircraft missions. 
Sea Ice Radar Experiments (SIRE) In 1978 and 1979, designed 
specifically to Investigate the effect of sea Ice ridging on radar 
Instruments. SIRE 1978 and SIRE 1979 occurred during the winter Ice 
seasons of those years. Both missions were conducted In the Beaufort 
Sea off northern Alaska, and the data for this thesir were collected 
exclusively In the Prudhoe Bay region. Typical aircraft ground tracks 
are shown In Figure 2.1. These flight lines consisted of flying over 
heavily ridged Ice at an altitude of 300 to 300 meters, for distances 
of approximately 25 km. During each flight line, data were obtained 
simultaneously from optical and radar sensors. 

This chapter will discuss only those Instruments used during 
SIRE which are of Interest to our Investigations; namely, the aerial 


BEAUFORT SEA 



Figure 2.1. Flight Lines for SIRE 1979 for Harch 12 
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mapping camera* the laser profiler and the microwave scatterometer* 

The first two are used to establish the Input data set to the radar 
backscatter model; whereas the data from the scatterometer are used 
for comparison and evaluation. Figure 2.2 Illustrates the positions 
of the Instruments on the airplane and the locations and relative 
sizes of their footprints (area of Illumination on the surface). Also 
discussed In this chapter are the data processing methods used to 
prepare the input to the model. Examples are presented which Indicate 
how the data was edited to constrain the ice surface parameters. 

2.2 Mapping Camera 

If one eliminates observations from ships* which are a form of 
"remote sensing*" then aerial photogrammetry Is the first remote 
sensing method used In sea Ice study* and certainly the first to 
provide wide area coverage. Although use of the .-•erlal mapping camera 
Is limited to clear weather and daylight* the Information available 
from this Instrument Is extremely helpful to Ice Investigators. 

Through photogrammetrlc interpretation It Is possible to 
calculate ridge height from shadows and to estimate age and thickness 
of the ice* provided that the snow cover has not obscured the surface 
texture and color tone. Figures 2.3 and 2.4 are facsimiles of typical 
aerial photograms. It is seen from these that the thinner the Ice* 
the greyer It appears In the photogram* with water and extremely young 
Ice appearing almost black. Multl'*year Ice (Ice which has survived at 
least one melt season) is much thicker than first-year Ice and can be 
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Figure 2.2. Relative Locations of Instruaients* Illtinlnatlon Areas 
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Figure 2.3. Facsimile of Aerial Photogram with Young Ice. Photo #04-181 
NASA JSC 375, courtesy NASA Langley Research Center 
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Figure 2.4. Facslnlle of Aerial Photogram with Multi-year Ice. 

Photo #10-031 NASA JSC, courtesy NASA Langley Research Center 
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differentiated from thick first year Ice because the former has a 
smoother, weathered and wind worn appearance* 

Because of its Increased accuracy, photogramme try Is used In 
various remote sensing experiments to provide the comparison data set 
(“surface truth**) against which Instrument responses are compared* 
Recent studies between data collected by surface based methods and 
photogramme trie data of the same areas have shown that even the 
experienced Ice observers In surfaced based camps consistently 
overestimated the Ice concentration and amount of ridging by 15-20%, 
and exaggerated the areas of older Ice by up to 40% [2]. 

The aerial mapping camera used during SIRE was the Zeiss RC-9 
with a 6*' focal length and a 74" field of view. This provided a 
ground coverage distance of 1*5 x the altitude (or 450 m at an 
altitude of 300 m) * Photograms were 9** square, black and white 
contact prints and were taken with 20-60% overlap to allow stereo 
analysis of the surface topography* First, the photograms taken with 
this camera were used to classify the Ice by age and thereby eliminate 
all non-first -year Ice from the data set; and secondly, they were used 
to verify that both the laser and the scatterometer were looking at 
the same features, even though their fields of view did not coincide* 
Also, stereo analysis was performed on a limited set of photograms to 
investigate the error Introduced In estimating a two-dimensional Ice 
ridging process using a one -dimensional model* 
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2.3 Laaer Profiler 

Another rcaote sensor^ the later profiler, hat been uted 
extenalvely froa aircraft to provide quantitative data on Ice 
roughnett. The profiler, utlng a radio frequency aodulated light 
beaa, aeaturea the phate delay between the tranaaltted and reflected 
beaa and thereby deteralnet the dittance froa the aircraft to the 
surface. Thut, when the aircraft filet at conttant altitude (relative 
to a«.an tea level), the dittance between the aircraft and the Ice 
surface Is a measure of the topography. Like the camera, the later 
cannot penetrate clouds or fog and It therefore limited In its use for 
general Ice monitoring. However, It has served Ice scientists well In 
the past by providing the Ice topography over extended areas. 

The profiler used In the SIRE missions was a Spectra Physics 
Geodollte Model 3A, (Figure 2.5), which contained a helium-neon laser 
(6328 A*) and a 8" telescope [lA]. The Instantaneous nadlr-looklng 
(0° Incidence) laser spot diameter was 7.5 cm at 300 m altitude; but 
the effective along-track resolution was degraded to 1.4 m (for 
typical velocity) due to the distance traveled by the aircraft during 
the Instrument's 20 ms Integration period. The relative distance 
measurement accuracy was typically 3-10 cm for smooth Ice [15]. 
However, errors could be significantly greater over ridges since 
ridges are comprised of piles of Individual blocks of Ice which can 
cause discrete distance steps rather than a continuous change In 
surface elevation. An analog output voltage proportional to relative 
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Figure 2.5. Laser Profiler, shown mounted Ir C 4 Ui»era bay of C-130, 
courtesy NASA Langley Research Center 
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height (typicelly 40 foot full ecele) wee digiclced et 100 He; and 
after soae preproceecing, which will be dlecuesed later In the 
chapter » these heights were used as the sole Input to the model* 

2 . 4 Microwave Scatteroaeter 

A microwave scatteroaeter Is an active (RADAR) Instrument which 
measures the absolute backscattered power from a target* When this 
target Is a surface of uniform scattering statistics, the use of the 
radar equation (see Appendix A) permits calculation of the surface 
reflectivity expressed as the normalized radar cross section, CJ* . 

The scatteroaeter was first used over Ice In 1969 by Rouse 
[16], who determined that It could be used to categorize sea ice based 
on age* The early success encouraged continued investigation, and It 
is presently known Chat 0** is a function of the electromagnetic 
parameters frequency, polarization ancl incidence angle, and of Che Ice 
parameters salinity, temperature and roughness [17,18]* 

The scatteroaeter used In the SIRE missions was the NASA 
Advanced Applications Flight Experiment Radiometer /Scatteroaeter (AAFE 
RADSCAT) [19], shown mounted in the aircraft in Figure 2*6* It 
operated at 13*9 GHz, transmitted a horizontally polarized signal, and 
received horizontally polarized backscatter at cn Incidence angle of 
48** The RADSCAT had a "pencil beam" antenna of beaawldch 1*25*, 
which, for altitudes of 441 m and 48* incidence angle, provided an 
instantaneous footprint on ttui surface of approximately IS m x 22 m* 
The RADSCAT used an "integrate and dump" method of data collection. 
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integrating for 100 ms, which corresponded to the Nyqulst spatial 
sampling for the antenna footprint* This Integration yielded a 
"smeared" footprint of 15 m x 29 m> Specifics on the Instrument 
characteristics for the RADSCAT can be found in Appendix A* 

The data from the scatterometer are used as the comparison data 
set for correlation with the model's predicted power* Once again, if 
the two data sets are highly correlated, then we can conclude that the 
scatterometer Is a useful sensor for the identification of Ice ridging* 

2.5 DATA 

This section discusses the differences between the Instrument 
perspectives of the Ice due to their differing nadir angles, and 
Identifies how these differences will be handled* It defines how the 
data set was determined and discusses the processing of the laser data 
which was necessary prior to using It In the model* 

2.5*1 Differences Between Instrument Perspectives 

As can be seen In Figure 2*2, the nadlr-looklng profiler and 
the aft-looking scatterometer have different perspectives of the Ice* 

In addition to their different spatial resolutions, the sampling rates 
are not the same, and their footprints are of different magnitudes* 

The laser's Integrated spot size Is 0*07 m x 1*4 m, while the 
RADSCAT' s Is 15 m X 29 m* Thus, the RADSCAT measures scattering from 
a two-dimensional area but the laser, because of the small width of 
Its spot, effectively sees a one-dlmenslonal view* The effects of 
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these differences must be compensated for In the model or be 
considered In Interpreting results* 

The difference In the viewing angle can be eliminated by 
rotating the coordinate system of the laser. It can be shown that 
this can be accomplished by adding the Incidence angle of the 
scatterometer to the surface slope as determined from the laser. The 
* difference In the sampling rates and spot length Is accounted for In 
the model by a weighted averaging of the laser data to approximate the 
scatterometer' s footprint length (alongtrack) . The consequences of 
using a one -dimensional representation of a two-dimensional scattering 
process will be discussed with the conclusions In Chapter Five * 

2.5.2 Data Selection 

The data set was established by locating on the photograms the 
footprints for the laser and the RADSCAT, and then making subjective 
Judgments that both sensors were "seeing” the same feature. The times 
for the areas used In this analysis were approximated from the 
photograms and Identified exactly from the lime series records of each 
sensor. Minor time adjustments were made between the laser and the 
scatterometer by matching dominant features In each time series. This 
was necessary for two reasons: each sensor was recorded using ' 
separate time reference and» more Importantly, the separation distance 
(time) between the sensor footprints was variable along the flight 
line because of changes In the aircraft attitude (roll, pitch, and 
yaw), ground speed and altitude* 
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For our analysis, all experloental measurements have been 
edited to select those data for which the degree of ridging Is the 
only significant parameter. All areas studied are strictly first-year 
Ice as Identified by photo-lnterpretatlon; they are, therefore, 
similar In salinity and thickness. Data from different years were 
obtained In the same general location and under similar ambient 
conditions as measured by on-board sensors and are believed to be as 
Identical as data from first year Ice along different flight lines In 
the same mission. 

Further, while It is possible that older first year ridges may 
be slightly less saline than the surrounding flat Ice^ (because of 
leeching and ablation). It Is not believed that this minor change In 
the Ice dielectric properties Is sufficient to affect the measure 
of significantly [20]. Only areas determined to be relatively free 
of loose snow were Included In the data set. The electromagnetic 
parameters affecting (f were held constant by eliminating any data 
from the study set which were not at horizontal polarization or 48* 
Incidence. 

2.5.3 Surface Slope Processing 

The Input to the radar backscatter model Is the mean Ice 
surface slope over distances which are small compared to the spatial 
resolution of the scatterometer measurement. This mean slope was 
obtained from the laser profiler's signal which contained Information 
on the Ice topography. Unfortunately, the raw laser data also 
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contained Instrument related artifacts (phase shifts) and distance 
changes due to aircraft vertical motion. Both of these perturbations 
had to be removed from the data prior to use In the model. 

A phase shift occurs every time the relative distance goes 
through a full scale change (typically 40 foot scale). Figure 2.7 
Illustrates how the laser signal must change phase to keep the signal 
on Its 40 foot scale as the aircraft changes altitude. 

A phase shift was removed by adjusting all subsequent data by 
the amount of the shift, which is dependent on the scale setting and 
the direction and magnitude of the shift. Shifts can be 25Z, 50%, or 
100% of the scale setting. For Instance, a 25% phase shift on the 40 
foot scale would require a correction of ten feet. 

Phase shifts are detected automatically by a program which 
computes the difference In magnitude between adjacent data points. 
Whenever this difference Indicates that a shift has occurred, the 
subsequent data Is Increased or decreased by the precalculated 
amount. Since the Ice surface does not change Instantaneously by ten 
or more feet, there Is little danger of editing out ridges with this 
program. Any jumps which are not detected by the program are removed 
by hand input corrections. Figure 2.8 shows the phase shift removed. 
It Is obvious from the figure that missed phase shifts would be 
quickly noticed In the data. 

Aircraft altitude changes are reflected in the laser data since 
the laser measures the distance between the aircraft and the surface. 


LASER STRIP CHART 



Figure 2.7. Physical Representation of Phase Shift 
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Figure 2*8. Removal of Phase Shift from Data 


Figure 2.9 llluecratee the aircraft motioa in the data. These slowly 
changing oscillations in the data could ideally be removed by a 
high-pass filter. Unfortunately, the spectra of the altitude change 
due to aircraft motion and of the ice topography overlap [21]; 
therefore high-pass filtering would also remove low frequency ice 
information, which would be unsatisfactory. Figure 2.10 shows the 
results of a high-pass filter attempt. Notice that the ridge sits in 
an artificial depression which affects the accuracy of the height 
measurement . 

The process which is used for this editing is a modified 
high-pass filter and substitution process which provides more accurate 
ridge information. The laser data (Figure 2.11a) is first low-pass 
filtered to provide the first approximation of the aircraft motion, 
(Figure 2.11b). The filtered line is then subtracted from the 
original data and the magnitude of this difference is monitored. 
Whenever the difference is greater than one meter, a substitution is 
made for the original data. This modified line now has had all 
significant ridging removed and replaced with flatter segments (Figure 
2.11c). The modified line is low-pass filtered once more to yield an 
improved approximation of the aircraft motion (Figure 2. lid). This 
second approximation is then subtracted from the original laser line 
to produce a flat and level sea level base reference for the surface 
profile (Figure 2. lie), which is the way we believe the ice freezes* 
The accuracy of this method could possibly be Improved by a second 
iteration of comparison between the first and second approximations. 



Figure 2.9. Aircraft Hotlon In Laaer Data 
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Figure 2.11. Aircraft Motion 
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Howtv«r» OM substitution hss providsd sn error of less then 5Z of the 
ridge height which was felt to be within accepteble limits for our 
study* it should be pointed out that the slopes to be calculated 
later use a differential height between points and this error in 
relative height has even less then effect on the slope* 

At this point the laser data is ready for input to the model* 
The next chapter discusses the modeling process* 




Figure 3.1. Onatott^s Signature Curvea 





33 


rad«r croat •action, 9*. ■■ a function of tha antanna angla, 6 (tha 

n 

•lavftlon angla of tha antanna aa Maaurad from nadir) •* Slnca flan 
lea aurfacaa axhlblt highly apecular (non-dlffuaa) acattarlng, tha 
•trongaat ratum occurs whan tha nadir angla la aaro dagraaa (0^), and 
tha atrangth of tha ratum signal dacraaaas aa tha nadir angla 
Incraasaa, aa llluntratad In Figura 3.2. Nota also fron this figure, 
that for flat lea, ^ • whara ^ la tha true Incidence angla and 

la defined aa the angla between tha direction of propagation for tha 
Incident wave and the local nonaal to the surface. Figure 3.3 shows 
chat for sloped surfaces, d > general, given the slope 

angle, a* of the surface, 

d ■ + a (3.1) 

If Ice ridges ware simply saooth, tilted plates of ice, then 
their reflective charactarlstlca could be deteralned froa the sig- 
nature for flat lea. Howavar, ridges ara not smooth, but are composed 
of blocks of differing slzas heaped together as the ice deforms. 

Figure 3.4 shows facslmllft of photographs of ridges taken on tha lea 
surface, and Figure 3.5 Is a facsimile of an aerial phoCograa of rough 
lea. Both figures Illustrate the chaotic arrangement ot the Ice 
blocks. Given this roughness which occurs In ridging, and all the 
Individual scattari.ug faces, tha assumption of specular scatter for 
ridged ice s&ema Inappropriate. Thus, a different signature curve 
(shown lu Figure 3.i) Is used for our model for ridged areas. Aa 
before, the cross section, a*, decreascc with iucreaslng nadir angle, 
but hare Che slope Is much shallower than that for flat lea. 
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Figure 3.3. Illustration of Equal Return for 
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Ridge height, 1.5 m, block thickness, 0.1 m 



Ridge height, 2.0 m, block thickness, 0.4 m 


Figure 3.4. Fecslmlles of Surface Photographs of Ridges, 
courtesy US Army CRREL 





ORIGIWi. y,-a^ 

OF POOR QUALITY 


Figure 3.5. Facsimile of Aerial Photogram of Ridged Ice, Photo #05-370 
NASA JSC 396, courtesy NASA Langley Research Center 
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3.3 Backacatter Model 

We see froe the signature curves In Figure 3.1 that the 
strength of radar backscatter. , is a function of tl^ ice surface 
roughness and the nadir angle of the antenna. The other surface 
parameters which are known to affect such as age, temperature and 
salinity, are held constant by editing the data selected for this 
investigation, as discussed in Chapter Two. Also, since the nadir 
angle for the RADSCAT measurements during SIRE was constant at 48° , 
the radar cross section for these regions of first-year sea ice should 
be a function of surface roughness only. Therefore, it is this 
information which is used as the input to the model to calculate the 
radar return. 

For this thesis the surface roughness Is provided by the laser 
profiler. The radar backscatter model uses this oae~dlmenslonal 
representation of the ice topography and Interprets it as subregions 
of slope (facets). Figure 3.6 Illustrates how one scatterometer field 
of view Incorporates many of these subregions, each of which lias 
associated with it an incidence angle, <^, and a strength of return 
power, or[(^]» as determined from the signature curves. The total 
return power for a footprint can be determined by summing the partial 
returns from each of the facets. 

The model is developed via several processes as shown in the 
flow chart in Figure 3.7. These processes can be broadly separated 
into the following: calculation of the incidence angle of each facet, 
<^(i): calculation of the return power from each facet, a[<^(l)]; 


HiaiMwvag 



Figure 3.6. Illustration of Facets In Equivalent Resolution Cell 
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Igure 3.7, Flowchart of Modeling Processes 
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spatial and teaporal averaging; and the addition of white noise 
(gaussian) • 

3.3.1 Calculation of Incidence Angle 

Radar return, or backacatter, is a function of surface 
topography. The laser profiler data, after the preprocessing 
described in Chapter Two to remove aircraft altitude perturbations, 
provides this input to our model. The profiler produces a time series 
of surface height, h(l), above a mean sea level reference, as shown in 
Figure 2. lie. Knowledge of the aircraft velocity (148 knots) and the 
sample rate of the profiler (100 Hz) allows calculation of the 
horizontal distance, d, between samples. With these values, h(l), and 
d, the slope, a, of each facet can be easily determined by the equation 

a(i) ■ arctan [h(i) ~ h(i-l)j (3.2) 

d 

In this way we have converted a sequence of height data in meters to a 
sequence of slope data* 

The slopes, a(l), are not equal to the incidence angles, (})(i), 
since the scatterometer and the laser did not have the same nadir 
angle. Equation 3.1 can be used however, and for 8^ ■ 48" we find 

<Kl) - a(i) + 48“ (3.3) 

This addition yields the same perspective of the ice surface that the 
RADSCAT had since it can be shown to be equivalent to the rotation of 
the laser's coordinate system. Figure 3.8 Illustrates, for a 
simplistic view of a ridge, the incidence angles for a back face. 



Figure 3.8. Illustration of 
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front f«c« and flat lea aa seen froa 48* nadir angle* After the 
addition in Equation 3*3 we have a time aeries of true incidence 
angles* defined with resptet to the RADSCAT's perspective, and still 
at 100 Hz. 

3.3.2 Instantaneous Reflected Power Model 

At this point we have an incidence angle, 4>(i), associated 

with the 1^^ facet in the field of view, and we use Onstott's 

signature curves [22] (Figure 3.1) as the basis for determining the 

return power, a[(^(l)], associated with that angle. The a{^(l)j are 

modified slightly from Onstott's in that they are a function of the 

incidence angle and not of the nadir angle* Secondly, they are linear 

rather than in dB since, as Indicated in the flowchart, the values 

will be averaged. Neither modification alters the way the curves are 

used, but merely simplifies the calculations* 

If the laser profiler had the same sample rate as the RADSCAT, 

as well as the same field of view, it would now be possible to compare 

the calculated returned power, Pgj^g|(l) (which at this point is 

proportional to <7[<^(1)]), with the measured RADSCAT value, Pg/g(l)* 

However, neither is tlM same for both Instruments. 

To compensate for the difference in the fields of view we can 

average over the length of an Instantaneous RADSCAT footprint all the 

9[^(l)j within that footprint. That is, since the instantaneous field 

2 

of view of the RADSCAT is 22 m (altitude x beamwidth / cos 6 ) and 

n 

the facets are 0.73 m (velocity / sample rate) apa.rt, there are 30 
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fac«cs within a glvan Instantanaoua footprint, aa shown aarllar In 
Figure 3.6* The reflactad power froa the surface illualnated by that 

Instantaneous footprint could now be represented by * 

30 

Pgia - 1 X a[(^(D] (3.4) 

30 i.i 

Recall froa Chapter TWo that the RADSCAT does not yield 
instantaneous reflected power, but rather Its output Is the result of 
an ’‘Integrate and dump” measureaent process. Thus, to calculate a 
returned power that Is coaparable to p^^^^ the averaging In Equation 
3.4 is not sufficient. It is also necessary to average over the same 
integration period that the RADSCAT does. This takes us to the third 
process in the aodel. 

3.3.3 Average Reflected Power Model 

Since the aircraft Is moving during the 100 as that the RADSCAT 
Integrates, it is not possible to deal with an Instantaneous field of 
view, as was done In Equation 3.4. Rather, we must force the model to 
average all the facets which are In the smeared (Integrated) area. 
Thus, we would expect the Index In Equation 3.4 to Increase, which It 
does. However, another factor must be considered — the radai antenna 
pattern. Like most antennae, the RADSCAT' s does not have a uniform 
response across the footprint, so we perform a weighting of the 
facets, to simulate this pattern across the smeared footprint. One 
last manipulation of the data takes place before the power Is 
averaged, and Is done for the sake of simplifying the required 
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calculations. This last procaas dlvidas tha ranga of poaslbla ^1) 
values into 4* blna and calculates tlw percentage of footprint in each 
bin. Once again, a flovchart. shown as Figure 3.9, is used to 
organize these processes. They are: the equivalent resolution cell, 
antenna weighting, percentage area for each incidence range, and 
summing the power. 

3. 3. 3.1 Equivalent Resolution Cell 

In the 100 ms the RADSCAT Integrates, the antenai has typically 
moved 7.3 n (velocity x integration time). This Increases the 
illumination area to 29.3 m in length, or the equivalent of 40 
facets. For all practical purposes, we can consider the time it takes 
the RAOSCAT to clear the Integrator to be negligible, and the next 
footprint area is immediately Integrated. Figure 3.10 shows three 
consecutive footprints to Illustrate the overlap which occurs on the 
ground. To simulate this, the model uses a "sliding** c.ata window 
wtiich processes 40 facets, slides the 10 facets equivalent to the 
100 ms integration wlme, and then processes the next 40 facets. This 

could be expressed, assuming no antenna pattern effect, as 

20 

Pslm(J) - -ri- ^ a[(^(10(j + 1) + 1)] (3.5) 

^1—19 j-l,2,...,K 

The time series p , (j) is now at 10 Hz, the sample rate of the 
Slfll 


RADSCAT. 
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Figure 3.9. Flowchart of Averaging Proc 
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Figure 3.10. Illustretlon of Overlapping Feaolutlon Celia 
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3. 3. 3. 2 Antenna Weighting 

The main beam radiation pattern for the RADSCAT Is reproduced 

In Figure 3«11 [23] • Due to this pattern and the movement of the 

antenna during the Integration period, the return power Is more 

strongly Influenced by the mld'-reglon of the Illumination area* To 

determine the weighting necessary to simulate these Influences In the 

model, a convolution Is performed* The Integral which represents the 

reflected power Is, 

[ +B/2 2 

W(0) - I G^[ T J F(6 - T) dT, 

J-B/2 

where G(6) Is the gain pattern (squared In the Integral In agreement 
with the radar equation), F(0) is the surface reflectivity as seen 
through a window of length B/3, and B Is the 3dB beamwldth of the 
antenna* It is consistent with work previously established by Beck 
[24]. For the RADSCAT, the antenna pattern Is approximated by 
G(0) - 1 + (-1*6167)(0^) + (l.O8O3)(0^) ~ (O.2837)(0^) 
for |0{ < 2*, and F(6) Is constant based on the small angle 
approximation that the variation of surface reflectivity with angle Is 
negligible* The result of the graphical convolution Is shown In 
Figure 3*12 and the continuous function W(0) Is approximated by the 
discrete weighting sequence, w^, shown as a dotted line* The figure 
shows that a simple weighting Is sufficient for the RADSCAT and 
amounts to doubling the power for the middle 20 values In the 

equivalent resolution cell. The returned power Is now 

20 

Pgia(j) - ^ S a[<|)(l0(j + 1) + 1)1 Wi (3.6) 

1—19 j-1,2,... 

where Wj^ ■ 2 for the closed Interval 1 ■ [-9, 10], and w^ • 1 
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Figure 3.12. Illustration of W«>lghtlng Sequence 
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elsewhere. At this point the series could be compared with 

the measured sequence from the RADSCAT on a polnt^per-polnt comparison. 

3. 3. 3. 3 Calculation of Area 

To simplify the calculations and because the slope of the 
signature curves Is small, the range of (f)(1) Is divided Into 4** 
subranges, or bins. Flat Ice Is defined as Ice with less than |2* | 
slope, l.e., an Incidence angle between 46-S(f . All other (^ are 
defined as rough/rldged ice. There are k 4^ bins. 

To calculate the returned power for an equivalent resolution 
cell, we determine what percentage of the cell Including the results 
of the weighting Is In each Incidence angle subrange. Once these 
percentages, p(k), are calculated, the return power for the cell Is 
determined by summing the product of p(k) and a((f>^), where a(<f>j^) 

Is the strength of return from the signature curve taken as the 
average value of 0 over the k^^ 4® bin. That is, for one equivalent 
cell, the modified calculation becomes 

Psim “ P(*^) <7(4>k) (3.7) 

k 

3.3.4 Addition of Radar Measurement Noise 

The signature curves provide 0(<|>^), which Is the mean power 
for a surface. It Is necessary that this me' el be capable of 
predicting the Instantaneous value, since the RADSCAT measures the 
estimated mean power which has been corrupted by noise. If we want to 
predict the Instantaneous power at the receiving antenna, we would add 
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Rayleigh noise* which would account for the fading statistics from the 
rough surface* If we are Interested In the received power measured by 
the RADSCAT (at the output of the scat terome ter), we would add white 
noise (gausslan) since a Rayleigh distribution becomes gauaslan when 
It undergoes the square law detection and Integration process. Since 
the comparison which needs to be made Is between the calculated power 
and the measured power, the white noise which Is added to the time 
series should have a gausslan distribution* 

The noise added In Che model had a zero mean, ai»l a standard 
deviation of 0.317, as determined from the bandwidth and Integration 
time for RADSCAT* After the addition of this noise, we have a time 
series which can be compared point -by -point with the measured values 
from the scatterometer* This comparison and the discussion of the 
results of the modeling effort will be presented In Chapter Four* 
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CHAPTER FOUR 


MODEL SIMULATION 


4.1 Introduction 

This chapter presents simulation results of the backscatter 
model, and comparisons with RADSCAT measurements obtained during 
SIRE. The selected data was divided Into two groups: a training set 
and an evaluation set. The processes for the model, as discussed In 
the last chapter, were used to calculate the simulated return, 
and then the measured return, from RADSCAT was used to “train" 

the model. That Is, using an Iterative process, the various 
components of the model were adjusted using a least squares approach 
to minimize the difference between the calculated values and the 
observed backscatter \ralues. Next, the refined model was compared to 
Independent data In the evaluation set. 

Also presented are the results of the stereo analysis which 
give Insight Into the two -dimensional nature of the ridged Ice 
surface. The chapter begins with a review of the training and 


evaluation data sets selection. 
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4.2 Data Selection 

Since the RADSCAT, laser profiler and RC-9 mapping camera were 
collocated on the same aircraft. It Is possible to use the measured 
aircraft drift to locate the laser spot and RADSCAT antenna footprint 
on the aerial photograms. Drift, or yaw (the difference between the 
aircraft heading and course), was recorded by the on-board naviga- 
tional system. By examining the traces of the footprints on the 
photos, a subjective decision was made to determine whether the laser, 
which tracked along the center of the photos, and the entire scatter- 
ometer cell were within an area homogeneous In roughness, l.e., seeing 
the saDM feature. Figure 4.1 shows an example of this process. Time 
was approximated from the photograms and determined exactly from the 
laser profile by feature matching. A number of areas were Identified 
and selected for processing In this manner, giving us confidence that 
the degree of ridging was the same for both sensors. 

4.3 Adjustment to <7[(^(i)] 

The radar backscatter signature curves shown in Figure 3.1 were 
derived by Onstott et al. [22], from Ice observations at the same 
frequency and polarization as RADSCAT. The ordinate is <7° , the 
normalized radar cross section, measured In dB. To review the 
notation used, Is the predicted total returned power for a 

given equivalent resolution cell, Pj^^g Is the measured returned 
power from the RADSCAT for one resolution cell. Both are linear 
values of power. For clarity, whenever It is necessary to use 



biFFEREfJr* 


DIFFERBIT TERRAIN 










56 


logarithmic values* they will be subscripted with dB« That la* 

Palm dB " ^°*10 Pglm* 

The term a[0] la the strength of return factor* determined from 
the cr* curve as a function of the incidence angle* It is used to 
calculate foi* this thesis* since an absolute ^"^^3 ( 

for RADSCAT) was not available. However* It was apparent after the 
first output of the model was analyzed* that vas not of the 

same order of magnitude as Pj(^y 3 * range of Pg^g^ was 10-50 

milliwatts* while the range of Pj^^g was 0-750 milliwatts* To 
normalize the simulated power to the measured power, It was necessary 
to find the proper scaling factor* m* and to scale c[(^] by that 
factor. That Is* 

0 ((^] - m X a* (((►] 

To determine m* an equality was forced between o® and ^ for 
smooth first year Ice* The values for a[(^] were then adjusted by m 
through the relation 

a(4)] - lO^^^dB ® (4.1) 

where m ■ 10 logj^Q m* 

4*4 Time Adjustment by Features 

Because of the scatterometer geometry* that Is* the 48* nadir 
arigle* there Is a delay of approximately six seconds between the times 
that the laser and the scatterometer view a particular feature. The 
time offset* t* can be calculated with the equation* 
t ■ (altitude X tan 48®)/ velocity 


I 


j7 

Since, for given flight lines, the altitude and velocity were fixed, t 
la nearly constant, and the adjustment Is easily made by using 
different start times In the computer program for each of the data 
sets* 

Occasionally, a different type of time adjustment must also be 
made to correct for the aircraft drift when the ridges are not 
perpendicular to the flight track* When the aircraft drift Is greater 
than 1*5*, the laser spot does not lie within the scatterometer 
footprint, and the two sensors cross the ridges at different points* 
Figure 4*2 Illustrates the difference In the distance between ridges, 
for each sensor, when the ridges are not parallel to each other nor 
perpendicular to the flight track* Since the comparison between the 
simulated and the measured values will be done on a feature-by~f eature 
basis, a time adjustment must be performed to prevent a degraded 
correlation* This correction Is made by compressing or expanding, as 
necessary, the time frame between the ridge responses for the 
simulated values* This correction does not alter the magnitude of the 
response* 

4.5 Adjustment to Signature Slopes 

When the new p . time series, calculated using the scaled 
Sin 

a[0] and AttKiL the time adjustment for features, was compared to 
Pr/S, It was determined that there was a need for further 
refinements to the model* Figure 4.3 shows the plot of the two time 
series* It Is seen from the plot, that the simulated power is too low 
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Initial Predicted Return and RADSCAT aerlea 



In rldgnd ar«a«» nnd too largo In araaa of aaooth lea* In addition, 
thara ia no variation in tha siaulatad aariaa in aaooth araaa. 
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An amilyaia of hiatograaa of povar varaua for both aariaa, 
ahovad that tha lack of reaolution in aaooth araaa waa cauaad by uaing 
only oiM bin for all flat ica. Thia producad an artificial grouping 
which affectivaly ignorad tha alopa of tha o* curva. Incraaaing tha 
nuabar of blna by raduclng tha ranga width to 1* providad adequate 
reaolution. 

To adjuat tha aiaulatad reaponaa for a batter fit to the aaxlaa 
and ainlaa of tha RADSCAT, it ia naceaaary to sonaider tha alopa of 
tha aignature curves. Referring to Figure 3.1, note that tha curve 
for flat ice is very non-linear, but over the narrow range 46-5(f it 
can be approxiaated by a linear expression. Utilizing tha 
slopa-lntarcapt fora, tha strength of return can be represented by 

<^f [♦JdB ■ ^ (4‘2) 

where is tha slope and b^ is the Intercept. There is a siailar 
approxlaatlon for a 

r ds 

With these approxlaatlona and Equations 3.7 and 41, p . can 

Sin 

be calculated as 

V in Kar(«k)dB/10) + (arC<l>k)dB/10) ] 

PslB “ («.3) 

k k - 1,2 20 

where is, by definition, zero for flat ice, k " (9, 10, 11, 12), 

and is zero for all other k. Tha separation of and 

characteristics is reasonable since flat ice and rough ice are 
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To find ctM boot flc of th« tvo Mrloo in thn cralalog soc, an 
itaratlva procaaa was davalopad tfhlch adjuatad and 

hg until tha quantity 

2 IPalaU) - PR/s(J)]^ 

i 

vas a nlulaua* That la, flndlnt tha laait squara arror by adjusting 
thasa four varlablaa producad tha final M>del* Flgura 4«4 shovt plots 
of a portion of tha training sat tiaa sarlas along with tha 
corraspondlng lasar traca and photograas* 

To avaluata tha quality of tlM fit batwaan tha slnul.Atad aiWi 
obsarvad sarlas, a statistical program davalopad by Sarvlca Dava 
Corporation (25] was usad on tha NASA coaputar systaa, a CDC 6600. 

Tha corralatlon coafflclant for tha training sat was .68. Thare ware 
716 data points, at 10 Hz, ylaldlng 90T confldanca bounda of 72Z and 
138Z. 

At this point, the "tuned" aodal was used on tha Independent 
evaluation set. All factors affecting tha model revalnad tha saaa. 

Tha only data aanlpulation on tha evaluation set was that for feature 
allgnaant. Tha correlation coefficient for the evaluation set was 
.63, using 729 points at 10 Hz. 

• 4.6 Stereo Analysis 

Recall that the profile provided by tha laser Is effectlvaly 
ona*^iaanslonal because of Its narrow call (.07 a), while the 
scatteroaetar has a resolution call which is approxlaataly 15 a vide. 
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The question which must be addressed Is, can one laser profile 
somewhere within the footprint's %rldth be truly representative of the 
average roughness seen by the scatterometer. To answer this, 
photographic stereo analysis was performed on two photograms to 
provide five parallel Ice topography profiles within the antenna 
footprint. Figure 4.5 Illustrates the location of the stereo profiles 
within the resolution cell. 

One region (A) evaluated was a perpendicular crossing of 
Isolated ridges composed of significantly high ridges with smooth Ice 
to either side of them. The other region (B) was of significant 
roughness with one prominent ridge and little smooth Ice. Figures 4.6 
and 4.7 are facsimiles of the original photograms analyzed. Figures 
4.8 and 4.9 are plots of the profiles along with the corresponding 
laser traces. 

The stereo profiles for each region were pairwise correlated to 
determine their similarities. Tnese correlation values are presented 
In Table 4.1, and are also shown plotted in Figure 4.10 as a function 
of the distance between the profiles. Note that the correlation drops 
to less than 50% within 28 feet, or approximately half the width of 
the RADSCAT footprint. To further describe the Ice surface the mean 
and standard deviation of the heights were determined for each 
profile. These statistics are also presented In Table 4-1. Note the 
similarity of the statistics for the profiles In each area. 
Unfortunately, the laser profiles for these areas could not be 
correlated with the stereo profiles due to distortion Introduced Into 




Facsimile of Aerial Photogram for Region A. Photo #05-370 
NASA JSC 396, courtesy NASA Langley Research Center 
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Figure 4.7. Facsimile of Aerial Photogram for Region B. Photo #05-405 
NACA JSC 396, courtesy NASA Langley Research Center 
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Stereo Profiles for Region B, with Corresponding laser 







Statistics Cross Correlation 

REGION A Mean a rtos xL28 xL14 xCL xR14 xR28 
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the photography which caused a serious alignment problem between the 
two series* (Ice scientists currently addressing the question of the 
accuracy of the laser have been successful In aligning 
stereo --photograph profiles with the laser by using much longer time 
records [15].) 

As a final evaluation, the stereo analysis profiles, and their 
average were used as input to the model* The simulated returns due to 
each are shown as Figures 4*11 and 4.12, along with the p^^^g for 
each area* Again, alignment problems due to distortion In the 
photography prevented us from correlating Pj^yg with these* 






CHAPTER FIVE 


CONCLUSIONS 


5.1 Remark » 

The goal of Chla thesis is to develop a model which uses 
one-dimensional surface topography to predict the backscattered power 
from ridged, first-year sea ice. Based on the empirical results 
presented in Chapter Four, we believe this goal has been achieved, 
within the limits of statistical uncertainty for geophysical data. As 
was stated In the Introduction, we know of no work which has addressed 
the radar response to ice ridges. Existing backscatter models are 
theoretical in nature and are based on small scale (relative to 
wavelength) roughness, lliere are no established results with which to 
compare this model. It stands on its own merits as the first step in 
developing the algorithms necessary to interpret radar returns irom 
ice. 

The fact that the training set had a relatively high 
correlation coeffeclent (.68) gives confidence that the model is well 
tuned. That the correlation is also re'<atively high (.62) for the 
evaluation set supports the conclusion that the model reasonably 


ttm 
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slmulactt ch« response of the scetteroaeter. Further, since the input 
to the Bodel is surface roughness it Is leesonable to conclude that 
the scatteroaeter does respond to roughness and can be used in the 
future to Indicate areas of ridging. 

The results of the stereo analysis shot; 4:hat, while the 
statistics for the profiles of a homogeneous (as defined subjectively 
from photogrAms) area may be similar, their point-by~polnt comparison 
is not necessarily good. This is reasonable given the rendomness 
which is inherent to ice roughness, particularly non'-rldged rough 
areas such as area B in the analysis. The decrease in similarity in 
the profiles across widths comparable to the scatteroaeter footprint, 
as illustrated in Figure 4.11, indicate that a onc'dimensional profile 
is not totally representative of the roughness within the footprint. 
Based on the correlation values though, it is an acceptable means of 
establishing the surface topography. 

It is unfortunate that the stereo profiles could not be 
correlated with the laser or the RAD3CAT series for the same areas, 
and thus provide statistical support for the model. However, the 
plots of these series shorn in Chapter Four do show a visual 
agreement, particularly between the model output for the profile 
average across the footprint, which would indicate a reasonably tuned 
model. 


5.2 Future Research 


If research like that reported in this thesis is to continue it 
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will be necessary to find either a better way to determine the ice 
topography than a one-dimensional laser trace » or, since the use of 
stereo analysis has severe limitations, a way to verify that the 
one-dimensional trace can be used to establish the true topography 
across the width of the footprint, as well as along Its length. 

Another future research consideration Is to develop the 
algorithms necessary to quantify the Ice roughness given the radar 
response, since this Is one of the needs of Ice researchers. It will 
first be necessary to determine the parameter, or set of parameters, 
which can definitively describe the roughness of Ice: something which 
may best be left to the physical scientists who have a more complete 
knowledge of the ice and Its make-up. The microwave scientists can 
then develop a means of interpreting the backscatter in these terms. 

Future research needs to consider also that sea Ice, 
unfortunately, does not consist of discrete areas of first-year Ice, 
homogeneous In roughness. For the large resolution cells necessary 
from satellite measui.'ements, a means of interpreting the signature 
from areas of mixed ages and degrees of roughness must be 
established. Additionally, since, for the scatterometer , the response 
to multi-year ice is very similar to the response to roughness in 
first-year Ice, It will be necessary to employ another sensor along 
with the scatterometer to provide accurate Information about the Ice 
age. At present, there Is no way, without a priori knowledge of the 
surface, to differentiate between these two types of responses. 

O 

Hopefully, the processing of the radar data to an absolute 0 will 
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provide the distinction needed* Until then> however, the passive 
radiometer would be a compatible choice since their responses 


complement each other* 
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APPENDIX A 


Referring to Figure A.l, the following are defined; 

0JJ » Nadir Angle of the Antenna 

B ■ 3dB Beamwldtii of the Antenna, in radians 

R ■ Slant Range ■ Altitude 

cos 6n 

W ■ Width of the Footprint ■ 2 R tan(B/2) 

- 2 R (B/2) (small angle approximation) 

- R B 


L ■ Length of Footprint 


A ■ Area of the Ellipse 


From the Radar Equation [18] 
P^G^X^a 


"r- 


(47T)V 


(projection to normal) ■ Width 

cos6 

n 

_ tt(R B)(R B) ^ irR^B^ 

4 cos6 4 cos6 

n n 


R B 

COS0 


Where, Pr * received power 

Pj ■ transmitted power 
6 ■ antenna gain 

X “ wavelength 
R * slant range 

a ~ radar cross section, the area intercepting that amount 
of power which, when scattered isotropically, produces 
an echo equnx to that received from the object 
0 ^ » the normalized radar cross section, or scattering 
coefficient. 


For 


Then, 


q 

A 


A, defined as above, is the area of the footprint 


R 


p^g^xVa 

(4ir)V 


P^G^X^O*R^B^ir _ P^G^X^o“B^ 


(4Tr)^R^4cos9 


256Tr^R^co8 0 


And, 


^R - G^X^O*B^ 


P 2 2 

256ir R cos6_ 
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Referring to Figure A. 2, 

'■t - =1 'x 


p, • « p„ 
In c X 


For Calibrate, v . ■ G. P. a, *nd 

cal 1 in 1 

or, 

(G^ If a loss, Gj^ Is a gain, but is not long term stable) 


V « G. G Pv CL 
cal 1 c X 1 


Then, 


. O, G, P- 
V , ■ 1 1 c T 

cal ; 


(A.l) 


For Operate, 
Then, 

» V ■ 

op 

From A.l, 


op 


G G P 
1 R R 


Pt “ 

«i ^-1 "c 


1 In 


and 


From A. 2, 


P ■ G P 
^in R R 






(A.2) 


^ a. G. G a. V G 
op _ lie ■ 1 op c 


=1 =R ''cal =T 


''cal ®E ‘'T 


From the Radar Equation, 


Solving for , 


R 


-2.. 2 _2 „o 

G A BP 

256 R^cos 6 
n 


256 01 . V G R^ cos9 
1 op c 1 

2 2 2 

\al S S ^ ® 


Note that G^ has been eliminated 








